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We study the magnetization of a Heisenberg spin ladder using exact diagonalization techniques, 
finding three distinct magnetic phases. We consider the results in relation to the experimental 
behaviour of the new copper compound GuziCsHi^N^iGli and deduce that the compound is well 
described by such a model with a ratio of 'chain' to 'rung' bond strengths (J/ J') of the order of 0.2, 
consistent with results from the magnetic susceptibility. The effects of temperature, spin impurities 
and additional diagonal bonds are presented and we give evidence that these diagonal bonds are 
indeed of a ferromagnetic nature. 

PACS numbers: 71.27.+a, 75.10.-b, 75.10.Jm 
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There has been considerable recent theoretical inter- 
est in coupled chain systems for a variety of reasons: 
Firstly the systems provide an interesting step from 
the relatively well understood one-dimensional behaviour 
towards two-dimensional systems (i.e. a dimensional 
crossover); A second reason for interest lies in the un- 
usual and exotic behaviour exhibited by spin ladder sys- 
tems, for example a spin gap Q and on doping, hole 
pairing and a finite superfiuid density A third, and 
indeed the dominant, motivation for this article lies in 
the increasing number of compounds which can be well 
described by considering the behaviour of strongly cor- 
related electrons confined to coupled chains. The com- 
pounds {V0 2 )P 2 7 [| and SrCu 2 3 [| may be de- 
scribed by ladder spin systems and recently doping has 
been achieved in Lai- x Sr x Cu02.5 [^)- In this article we 
shall concentrate on the magnetic behaviour of the new 
copper compound C-^C^i^A^^C^ || which in con- 
trast to the other examples, exhibits a spin gap which 
is relatively small, thereby allowing a study of the mag- 
netic effects with a relatively modest magnetic field. We 
will show that the magnetization of the material is well 
described by the Heisenberg model on a ladder system. 
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FIG. 1. The Cu-Cu super-exchange paths in the com- 
pound Cu2 {CsHiiN'i^Cli (taken from Chaboussant et at). 



Recent experimental work on Cu2{C§Hi2N2)2Cl4 : , in- 



cluding magnetization, susceptibility and spin resonance 
experiments, has been presented by Chaboussant et al [pi. 
The material is thought to consist of effectively isolated 
coupled chains as shown in figure [l]. Superexchange gives 
rise to a coupling along the chains (strength J) and an 
interchain coupling (strength J'); there is an additional 
diagonal coupling J c {,. We shall firstly neglect J c i, since 
its relative strength is believed to be small although we 
shall include it later in the paper. Using the suscepti- 
bility data, perturbation theory and a high temperature 
series expansion, Chaboussant et al have deduced a bond 
ratio J'/J ~ 5.5. 

The Hamiltonian we shall use to describe the com- 
pound is the Heisenberg model on a ladder (2 x L) sys- 
tem, defined by 

U = J' ^2 S i.! ' S ^ 2 + J ^2 ' S J+b/3 

3 0,3 

+ J2 9a^ B H a Sl p (1) 

0,3,<x,l 



where (3 (=1,2) labels the two legs of the ladder (oriented 
along the x-axis), j is a rung index (j=l,...,L) and J and 
J' are the bond strengths along and between the ladders 
respectively. The final term represents an applied field in 
the direction a; we simplify this term to <?/ib H z S z al- 
though we should note that anisotropic effects may have 
minor but observable effects. 

The behaviour of this Hamiltonian in zero applied field 
(H = 0) and at zero temperature is now relatively well 
understood Q and is perhaps best understood by first 
considering the limit J = 0. In this case, the ground 
state has total spin zero and is formed by creating a sin- 
glet bond on each rung; excitations require one of these 
singlet bonds to be broken to form a triplet at an energy 
cost J'. This gapped state persists with the introduction 
of interchain coupling J and the triplets can propogate 
and form a coherent band with dispersion J' + Jcos k. A 
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gap, and the associated dispersion, has been observed ex- 
perimentally in the compound we are considering 0] and 
it is believed that in zero applied field the energy spec- 
trum of the system consists of a total-spin-zero ground 
state, a gap to the first excited state (triplet) and at 
higher energies, (bands of) states with even larger values 
of total spin. 

In this paper then, we will consider the effects of an 
applied field on the system i.e. considering the magne- 
tization curve. In the first section, we calculate A4(H) 
at zero temperature, and by comparison with experiment 
we deduce the ratio J/ J' relevant to the compound. We 
then consider the effects of a finite temperature, deduce 
the relevance of random spin impurities, and calculate 
the effect of introducing the small diagonal interaction. 
Finally, we discuss other possible relevant factors. 

The technique we have used in this study is Lanczos ex- 
act diagonalization on 2 x 12 and 2 x 16 ladder geometries 
with periodic boundary conditions. We have considered 
the momentum of the states (k x = =f-m where m is an in- 
teger) and also the parity of the states under a reflection 
in the symmetry axis along the ladder (even (R x = 1) or 
odd (R x — —I)). Since the hamiltonian commutes with 
the component of total spin in the z-direction (S z ), we 
may consider the subsets of S z individually. For a spe- 
cific value of applied field, we consider the lowest energy 
state for each subset and can then easily apply the field 
dependent term of the hamiltonian oc —H.S Z . With in- 
creasing applied field, states with larger S z become more 
favourable and we obtain a 'staircase' of states in the 
magnetization curve until saturation. 

In figure || we show results from the 2x12 system for 
various values of the ratio J/ J' (and also the case J/ J' = 
0.2 as an example of the 2 x 16 system). In addition, 
we plot the experimental results for the lowest available 
temperature (0.42K). The magnetization is normalized in 
such a way that the saturated system has magnetization 
unity and the applied field is normalized such that the 
value at which magnetization becomes non-zero is unity. 
For the cases with J/ J 1 equal to 0.5 and 1.0, we have 
plotted the staircase structure resulting from the finite 
system; using the midpoints of the 'steps' we have drawn 
a smooth magnetization curve. For the case J/ J' = 0.2, 
the results from the two system sizes are almost identical 
indicating that the finite size effects are very small (this 
is also true for other J/J')- 

Immediately we notice that two critical fields may be 
defined: For an applied field H < H c \ the magnetization 
remains zero (in a singlet ground state); the field H c \ cor- 
responds to the singlet-triplet gap at zero applied field; 
Then, with increasing magnetic field H c \ < H < H C 2, 
the magnetization increases until it reaches its saturation 
value at H c i- At this point it is worthwhile mentioning 
the work of Affleck |?| concerning gapped, integer spin 
antiferromagnetic chains: In an axially symmetric situa- 
tion the ground state above H c \ may be considered as a 
condensate of low energy bosons; varying the field varies 
the chemical potential of the bosons and the boson num- 



ber corresponds to the magnetization. In the limit of zero 
boson density (H \— > H c \ from above) the magnetization 
is shown to behave as M. (H) oc y/H — H c i and there is a 
power law decay in the staggered magnetization orthogo- 
nal to the applied field ||: For H c i < H < H c2 the spins 
exhibit a canted spin structure with a uniform magnetic 
moment in the direction of the applied field. The square 
root singularity appears consistent with the theoretical 
behaviour in figure 0. 
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FIG. 2. Magnetization as a function of applied field. 
A'! (sat) is the saturated value of the magnetization and A 
is the singlet-triplet gap with zero applied field. Results are 
shown for J/J' equal to 0.2, 0.5 and 1.0. along with the ex- 
perimental data (bold line). 

We emphasise that the theoretical results with J/ J' = 
0.2 are extremely close to the experimental data, and this 
ratio is consistent with that deduced by Chaboussant et 
al by analysing susceptibility data. Notice however the 
rounding of the experimental data in the region of the 
critical fields H c \ and H C 2 and the stronger singular be- 
haviour in the theoretical results. The aim of the remain- 
der of this article is to discuss the origin of these effects 
and we firstly extend our results to take into account the 
finite temperature, specifically looking at the rounding 
behaviour close to H C 2- 

It is easy to calculate the complete spectrum of en- 
ergy levels (by considering each subset of S z separately) 
and consequently the thermodynamic quantities can be 
calculated. For a specific value of the applied field the 
magnetization is defined by 

KA E» A exp(-/?i^)^ 

M = — 5 — (2) 

£ nA exp (-/?££* ) 

where is the energy of the nth eigenvalue of the 
hamiltonian (equation Py) with a z-component of spin S z . 
Since we are only interested in the region of the magne- 
tization curve close to H c i-, we restrict the summation 
over S z to S z > 8 (for the 2 x 12 ladder); states with 
smaller S z contribute only minimally in this region (this 
has been checked) and also subsets with S z < 8 include 



2 



many more states and computational limitations become 
important. 

The resulting magnetization curve for various values 
of (3 is shown in figure [| we concentrate on J / J' = 0.2 
since this corresponds closely to the effective ratio in the 
compound. It is immediately obvious that the effect of 
temperature is to cause a rounding of the magnetization 
curve, much as observed experimentally. From experi- 
mental considerations, we would expect a temperature 
corresponding to (3 ~ 4^ ~ jjj^j§ ^31 (where the 13.2 
originates from the fact that we have normalised such 
that J' is unity, see reference ||). Therefore, whilst tem- 
perature does indeed effect the magnetization curve in 
the vicinity of the critical fields, there must be some other 
factor to explain the small discrepancy between the the- 
ory and experiment. 
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FIG. 3. Magnetization as a function of applied field for 
various values of temperature (1//3). The experimental data 
(0.42K i— > f3 ~ 30) is shown as a bold line. 

The next step is the inclusion of random spin impuri- 
ties, i.e. including a term in the hamiltonian to describe 
the interaction of the local spins with random magnetic 
fields. The additional term included in the hamiltonian 
has the form 

T-tirnp = ^ W O,P S 3,0 ( 3 ) 

where Wj t p is the impurity strength at the site j, (3 (as 
defined in the initial hamiltonian), chosen randomly be- 
tween —w/2 and w/2 and S| p is the z component of spin 
on that site. In order to conserve the reflection symme- 
try we choose the impurity weight of the two sites on 
a particular rung to be equal (wj,i = Wj.?)- The inclu- 
sion of random spin impurities breaks the translational 
symmetry of the system (k x is no longer a good quantum 
number) and we also note that it is necessary to average 
over several realizations of the disorder due to statistical 
fluctuations. 

In Figure ^ we show the effects of including random 
spin impurities with various weights (w). The dominant 



effect of the impurities appears to be an increase in the 
critical field H C 2 ; The effect on the rounding of the mag- 
netization curve does not however appear to be the origin 
of the small discrepancy in shape of A4(H) between the 
experimental and theoretical results. 
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FIG. 4. Magnetization as a function of applied field in- 
cluding random spin impurities. Various impurity weights 
are shown for a temperature corresponding to (5 = 30. The 
experimental data is shown as a bold line. 

As a final investigation, we examine the possible im- 
portance of a diagonal cross bond (J c b) as shown in figure 
|l|. To include this affect we add the following term to the 
hamiltonian 

T~tcb = Jcb ^ Sj,l ■ Sj + 1,2 (4) 

i 

As suggested experimentally H, we take J c b to be smaller 
than both J and J' . In figuresHa and b we show the mag- 
netization curve for various sets of parameters, the choice 
of which has been led by the desire to keep H C 2/H c i close 
to the experimentally deduced value (increasing J c b in- 
creases the effective interchain coupling and to keep this 
ratio constant, J must simultaneously be increased). The 
results were calculated for zero temperature on a 2 x 10 
system (it has been verified that finite size effects arc neg- 
ligible). As previously, the results are normalised such 
that J' is unity, saturated magnetization has value unity, 
and H c i is unity. In figure ||a we show the experimental 
data and the theoretical results corresponding to J = 0.2 
(Jcb = 0) as shown previously, and J=0.225 with the 
cross bond strength ranging from to 0.2. A similar plot 
(figure ||b) shows data corresponding to J = 0.18 and 
negative cross bond strength ranging from to -0.15. 

Several interesting features are apparent in the results. 
Firstly, introducing the cross bond interactions does not 
affect the overall behaviour of the magnetization curve. 
The major effects are to firstly change H c i (this effect 
is not apparent due to the normalization), and secondly 
to change the shape of the magnetization curve slightly. 
Positive J c b infact changes the shape away from the ex- 
perimental behaviour, seeming to increase the singular- 
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ity behaviour. Surprisingly, a negative J c & appears to 
shift the theoretical curve closer to the experimental be- 
haviour, with say J = 0.18 J c b = —0.1 being a reasonable 
parameter choice. The results seem therefore to suggest 
a ferromagnetic diagonal interaction, and further inves- 
tigations of the orbital behaviour in the compound are 
required to deduce if this is reasonable. 
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may arise due to the fact that in the experimental data 
provided, the field H a is not applied perpendicular to 
the plane of the chains and hence the magnetization M 
is not parallel to the applied field: there is a g— factor 
anisotropy. 

Another point we should make is that we find various 
choices of parameters which give reasonably good agree- 
ment with the experimental results. Fitting the mag- 
netization curve alone is not sufficient to allow the rel- 
ative strengths of the parameters to be fixed extremely 
accurately (the Heisenberg chain with both nearest and 
next-nearest neighbour interactions || also gives similar 
results). We can however make reasonable deductions 
about the effective interactions present in the compound. 
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FIG. 5. Magnetization as a function of applied field in- 
cluding diagonal cross bonds J c f, of various strengths. Figure 
5a corresponds to J = 0.225 and positive J c b whilst figure 5b 
corresponds to J = 0.18 and negative J c (,. The experimental 
data is shown as a bold line and the result corresponding to 
J = 0.2 J c f, = is also included for comparison. 

Summarising the results above, we find that the 
Heisenberg model on a ladder geometry describes well 
the magnetization of CuziCsHizNz^Cli, with a ratio 
J I J' ~ 0.2. Temperature causes a rounding of the mag- 
netization curves in the vicinity of the critical fields and 
random spin impurities tend to increase H c i- Diagonal 
bonds improve the theoretical/experimental agreement 
further and we suggest that the diagonal interactions ap- 
pear to be ferromagnetic. The discrepancy that still ex- 
ists however is perhaps not surprising since we are us- 
ing only a simple Heisenberg model. Some complications 
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